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We calculate the nuclear enhancement in Drell-Yan transverse momen- 
tum distribution in hadron-nucleus collisions. In terms of multiple scatter- 
ing picture, we compute the contributions from double scattering. Applying 



Oh! 

1^ ■ the generalized factorization theorem, we express the nuclear enhancement in 



terms of twist-4 nuclear parton correlation functions. We demonstrate that 
these twist-4 nuclear parton correlation functions are universal, and are the 
same as those used to explain the nuclear dependence in di-jet momentum 
imbalance and in direct photon production. Using the known information on 
the twist-4 parton correlation functions, we estimate the nuclear enhancement 
for Drell-Yan pair in large qt region. We also discuss the source of nuclear 
suppression in small qx region. 
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I. INTRODUCTION 



Recently, it is found in experiment that the transverse momentum distribution of Drell- 
Yan pair show anomalous nuclear dependence The observed nuclear enhancement in 
average squared transverse momentum, (gf^), grows approximately as A^^^ with A the atomic 
number of the nuclear target |T|,0]. 

Since the invariant mass of the observed Drell-Yan pairs, Q^, is large, it is necessary to 
understand the nuclear enhancement at partonic level. Any single scattering with a large 
is localized in space, and hence, does not introduce the A^^^ type nuclear dependence. The 
fact that the observed nuclear enhancement grows with nuclear size clearly indicate that the 
anomalous nuclear enhancement should be a result of multiple scattering of a parton inside 
the nucleus. 

In terms of factorization at higher twist, Luo, Qiu and Sterman (LQS) have developed a 
consistent treatment of multiple scattering at partonic level 0. LQS expressed the nuclear 
dependence of di-jet momentum imbalance in photo-nucleus collisions in terms of twist-4 
nuclear parton distributions. Using the Fermi Lab E683 data, LQS estimated the size of 
the relevant twist-4 parton distributions to be of the order of 0.05 — 0.1 GeV^ times typical 
twist-2 parton distributions [Q. 

The transverse momentum distribution for Drell-Yan process, da/dQ'^dq'^, is a typical 
two-scale problem in perturbative QCD. When is the same order as Q^, da/dQ'^dq^ 
becomes an one scale process, and the conventional perturbative QCD calculation works 
well. However, when is much less than Q^, do jdQ^ dq^ becomes perturbatively unstable, 
and is proportional to l/o'l- as — * 0. In addition, the conventionally calculated partonic 
parts develop two powers of large logarithm, log^(Q^/g|^), for every power of as, due to 
collinear and soft gluon radiation from incoming beam partons. It is therefore necessary to 
perform the resummation of these large logarithms in order to get a correct Drell-Yan q^ 
spectrum for the region when q^ <^ Q'^ 0. 

In terms of multiple scattering picture, nuclear dependence of Drell Yan q^ spectrum is 
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a result of multiple scattering between incoming beam parton and nuclear matter before the 
Drell-Yan lepton pair is produced. When and are both large, double scattering should 
be dominant. Method developed by LQS can be naturally applied for calculating nuclear 
dependence of Drell-Yan da / dQ^ dq^. However, when is large and is relatively small, 
calculating the nuclear dependence of Drell-Yan da/dQ"^ dq^ will face the same complication 
as calculating da/dQ"^ dq^ itself. As emphasized in last paragraph, when g|, <^ Q^, Drell-Yan 
da/dQ"^ dq^ is not perturbatively stable without resummation of large logarithms. Due to 
possible interference of multiple scattering and coUinear and soft radiation, the resummation 
of large logarithms may not be as straightforward 0. 

Because of two physical scales, and even at lowest order, there are two types of 
double scattering: soft-hard double scattering and double-hard double scattering. In soft- 
hard process, the parton from the beam first absorbs a soft gluon at the amplitude level, 
and then undergoes a hard scattering with another parton from the nucleus to produce the 
virtual photon 7* at and q^. In double-hard process, the parton from the beam first 
undergoes a hard scattering with a parton from the nucleus to produce a quark q (or an 
anti-quark) at the order of then the produced quark q (or anti-quark) annihilates with 
an anti-quark (or a quark) to produce the virtual photon at and q^. As explained later, 
the double-hard process resembles the classical double scattering picture, while the soft-hard 
process does not. These two types of processes have opposite sign at the amplitude level. 
When g-r — i> 0, the amplitude of soft-hard and double-hard scattering cancel each other. 
Both the soft-hard and double hard process contribute to nuclear enhancement, but their 
interference gives suppression effect. In large g-r region, these two processes require very 
different parton flux, and the interference term is less important. However, in smaller qx 
region, the interference term becomes more important and can reduce the enhancement, or 
even give suppression effect. 

In this paper, we calculate the nuclear dependence of Drell-Yan da/dQ'^ dq^ for large gr 
region, where fixed order perturbative calculation works well. Nuclear dependence in small 
qT region will be addressed in Ref. P|. This paper is organized as the following, in Sec. II, we 
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outline the general formalism and consideration for double scattering contribution; in Sec. Ill, 
we compute the soft-hard double scattering contributions; in Sec. IV, we derive the formula 
for double-hard scattering process; in Sec V, we discuss the relation between soft-hard and 
double-hard processes; in Sec. VI, we present the numerical result and conclusions. 



II. GENERAL FORMALISM 

To study the nuclear dependence of transverse momentum distribution, we consider the 
differential cross section for Drell-Yan process: 

h{p')+A{p)^l+r{q)+X , (1) 

where p' is the momentum for the incoming beam hadron, and p is the momentum per 
nucleon for the nuclear target. In Eq. (P, g is the four-momentum for the virtual photon, 
which decays into the lepton pair. In terms of contributions from multiple scattering, we 
expand the differential cross section as 



dQ'^dq^dy dQ'^dq^dy dQ'^dq^dy 



(2) 



The superscript "S" and "D" denote contributions from the "single scattering" and "dou- 
ble scattering", respectively, and "..." represents contributions from even higher multiple 
scattering. In Eq. (^, Q is the invariant mass of the lepton pair, qx is the total transverse 
momentum of the lepton pair, and y is its rapidity 

1, go + qz /o^ 

1/ = o In . (3) 

2 qo-Qz 

As discussed above, single scattering is localized, and therefore, da^^/dQ'^dq^dy does not 
contribute to anomalous nuclear dependence. It is proportional to A, modular nuclear 
shadowing effect. On the other hand, clear A^^^ dependence shown in the experimental 
data indicates that multiple scattering (higher than double scattering) is less important 
in proton-nucleus collisions. In this paper, we compute the double scattering contribution 
daj^y^/dQ'^dq^dy. 



For fix target experiments, often tlie variable xp = '2qz/ ^/s is used. We can easily convert 
our results dahA/ dQ'^dq^dy to dahA/dQ^dq^dxp by using the relation 

dQ^dq^dy = , = dQ^dq'idxp . (4) 



In order to compare with experimental data, we calculate the ratio of total differential cross 
section and single scattering contribution. 



dQ'^dq^dy / dQ'^dq^dy dQ'^dq^dy / dQ'^dq^dy 



For single scattering. Fig. [l| shows the lowest order Feynman diagrams that contribute 
to transverse momentum distribution of Drell-Yan pair. According to factorization theorem 
0, the single scattering contribution can be expressed as 

da^ 2al^a^ 1 r dx' dx 



"^a;;^. - "7 V ^^/"^"^'^ ^'^/^^^^ '^^'^^^ 



dQ^dq^dy 3Q 

x's + u-Q-^ y^ x's + u-Q^ )■ ^ ^ 

In Eq. (^, x' is the momentum fraction of parton "c" from the beam hadron h, and x is 
the momentum fraction of parton "a" from the the nuclear target. fc/h{x') and fa/A{x) are 
the parton distribution functions for the hadron and the nucleus respectively, s, t, u are the 
invariant variables defines as 

s = {p + p'f = 2p-p' ■ (7a) 
t = {p'-qf = -2p'-q + Q^- (7b) 
u = {p-qf = -2p ■ g + g2 . (7c) 

The averaged matrix element square |Mafe^^. p can be computed from the Feynman diagrams 
shown in Fig. |l] p|. For the annihilation diagram 



|M„--rP = 5e^2(-. + ,+^j ; (8a) 

for "Compton" diagram with the quark from the beam and the gluon from the target 



for "Compton" diagram with the gluon from the beam and the quark from the target 

,9 1 9 ( s -u 2QH\ , , 

|M,,^,.p = -e -2 + ^ Sr ■ (8c) 



In Eq. d^), the factor of couphng constants e^gl is aheady included in the over all factor in 
Eq. s,t, u are defined as: 

s = {xp + x'p')'^ = xx's ; (9a) 

i = (x'p' - qf = X\t - Q2) + Q2 . (g^^ 

u = [xp- qf = x{u - Q^) + . (9c) 

According to generalized factorization theorem , double scattering contribution to the 
cross section can be factorized into the following form: 

dQ^d^d^ = ^ I dx'f,/h{x') , (10) 

where da^j^ represents the double scattering contribution between a beam parton of flavor 
"c" and the nuclear target. The fc/h{x') is the normal twist-2 parton distribution of flavor 
"c" within a hadron "h" . We can simplify the calculation by integrating over the leptonic 
part of phase space first. 

daj^A^a _ 1 e2 d&aA^r ^^^^ 



dQ^dq^dy 27r 3Q2 dQ^dq^dy ' 

According to the factorization at higher twist, da^^/ dQ'^dq^dy can be expressed in the 
following factorized form: 



d^cA 



dQ'^dq^dy 



J dxdxidx2Ti{x,xi,X2,p)Hi{x,xi,X2,p,p',q). (12) 



The graphical representation of Eq. ( P^ ) is shown in Fig. |^. 

At lowest order, there are two types of double scattering. Fig. ^ shows sample diagrams 
at amplitude level for these two types of double scattering. In double scattering process, 

6 



there are two partons "a" and "b" from the nucleus participate in the scattering process. The 
kinematics can only fix one parton momentum. We need to integrate over the other parton's 
momentum. As shown in Fig. because of the extra scattering, we have two propagators, 
that is, two possible poles. Therefore, the leading contribution of the integration over the 
extra parton momentum is given by the residues at these poles. At fixed Q and g-r, "we cannot 
have both of these two poles at the same time. The total double scattering contribution is 
the sum of the residue contributions from these two different poles. Correspondingly, there 
are two types of double scattering. Taking the residue corresponding to the pole of the first 
propagator, as shown in Fig. is equivalent to say that the first parton (gluon) is soft. We 
call it soft-hard double scattering. In this kind of processes, the parton "c" from the beam 
first absorbs a soft gluon at the amplitude level, and then undergoes a hard scattering with 
another parton from the nucleus to produce the virtual photon 7*. The large momentum 
transfer occurs only in the hard scattering. Absorption of the gluon modifies the effective 
parton fiux, so as the overall cross section. Taking the residue from the pole of the second 
propagator, as shown in Fig. correspond to a situation when both partons a and h have 
non-zero momenta. We call this kind of processes double-hard double scattering. In this 
type of processes, the parton "c" first have a hard scattering with a parton "b" from the 
nucleus to produce a quark q (or an anti-quark) at a scale of then the produced quark 
q annihilate with an anti-quark (or a quark) "a" to produce the virtual photon at a scale 
of Q^. More discussions about the difference and the relations of soft-hard and double-hard 
processes is given in Sec. V. In the following two sections, we derive the contributions from 
soft-hard and double-hard processes respectively. 



III. SOFT-HARD DOUBLE SCATTERING CONTRIBUTION 

In this section, we derive the contribution from soft-hard scattering. At lowest order, 
the soft-hard double scattering has three types of subprocesses, as show in Fig. ^, Fig. |^ 
and Fig. |^. These three types of subprocesses correspond to adding two soft gluons to the 
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lowest order "Annihilation" and "Compton" subprocesses. In the following subsections, we 
present the derivations for one subprocess, and provide the results for other subprocesses. 



A. Factorization and expansion for soft-hard processes 

Consider subprocess shown in Fig. ^. There are three independent four-momentum 
linking the partonic part and corresponding two-quark- two-gluon matrix element. In the 
center of mass frame of high energy collision, all partons inside the nucleus are moving almost 
parallel to each other, along the momentum of the nucleus. All three parton momenta can 
be approximately replaced by the components coUinear to the hadron momentum. After 
such coUinear expansion, the double scattering contribution for contribution from process 
shown in Fig. ^ can be written as: 

da^^/J , 1 f f - 



j dxdxidx2 J d'^kTT{x, xi, X2, kT,p) H{x'p' , X, xi, X2, kT,p, q) , (13) 



dQ'^dqj.dy 2x's 

where 2x's is the flux factor between the incoming beam quark and the nucleus, and x'p' is 
the momentum carried by the beam quark. We keep the /c^ for the soft gluons in order to 
extract a double scattering contribution beyond the leading twist. The superscript "S'if" 
denotes "soft-hard" double scattering. In Eq. (|13D, the two-quark-two-gluon matrix element, 
T, is defined as 

dy- dy^ dy^ d'^yr 



T(x,Xi,X2, kT,p) = J 



2n 2tt 2ti {2tiY 

X ^ {pA\A-^iy2, Ot) ^,(0) 7+ My) A^iVi^ Vit) \Pa) . (14) 

The corresponding partonic part H is given by the diagrams shown in Fig. with gluon 
lines contracted with p^p"^ ^ quark lines from the target traced with {'-f ■p)/2, and quark lines 
from the beam traced with (7 ■ {x'p'))/2. Here, we work in Feynman gauge, in which the 
leading contribution from the gluon field operators is ^ A'^{p'' /p'^). 

In order to extract the lowest order high-twist contribution, we expand the partonic part 
H introduced in Eq. (|^) at kx = 0, 
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H{x'p',x,xi,X2,kT,p,q) = H{x'p',x,xi,X2,kT = 0,p, g) 



+ 



OH 



1 d^H 



2 dk^dk^ 



Ua 7./3 _|_ 



(15) 



kT=0 



kT=0 lyrvrpurvrp 

In the right-hand-side of Eq. (p!5D, the first term is the leading twist eikonal contribution, 
which does not correspond to physical double scattering, but simply makes the single- 
scattering matrix element gauge invariant. The second term vanishes after integrating over 
/ct- The third term will give the finite contribution to the double scattering process. Sub- 
stituting Eq. ( [T3| ) into Eq. (|T3|), and integrating over (PkT, we obtain 

I 



dQ^dq^dy 2x's J 2n 2n 2n 2 



dy- dy-, dy, 1 ^^(q) ^+ F+"(yr)bA) 



X 



-2' 



a/3 



2 dk^dk^ 



H{y ,yi,y2,kT = 0,p,q) 



(16) 



In Eq. ([iq) , the modified partonic part H is defined as 



H{y-,y^,y2,kT,p,q) = J dx dxidx2e'''P''y e'^'^v+v, e-'^^i-''^)p+y2 

xH{x'p',x,Xi,X2,kT,p,q) , (17) 



where the partonic part H is defined by Eq. (|T^), and is given by diagrams shown in Fig. ^. 
In Eq. (|16]), = F^^np, and F^"" is the field strength, and vector rip = 5^+. In obtaining 
Eq. (P^), we use the factor k^k!^ in Eq. ( [T3| ) to convert the field A'^ into field strength by 
partial integration. 



B. Integration over the parton momentum fractions 

From Eqs. ([I6D and (p!?!), all integrals of x, xi and X2 can now be done explicitly without 
knowing the details of the multi-parton matrix elements. 

Consider the diagram shown in Fig. The final state photon-gluon two particle phase 
space can be written as 

1 1 , / i — k^ — 2kT ■ q\ , . 

Idtt^ x's + u — \ x's + u — Q'' J 
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In deriving Eq. (0), we have omitted the factor dQ'^dq^dy, due to the definition of the 
invariant cross section, = —kTak^, and t is defined in Eq. 

Using Eq. (|T8|), the contribution to H from the diagram shown in Fig. ^ can be expressed 



as 



1 - 

Hla = -^Ci — ■ — Hia{x, Xi, X2) 

4tt x's + u — 



1 i-kl-2kT-q\ 

Hx + xi + ^^^ i-^ , (19) 



xi — X2 — — it xi — + ie \ X s -\- u — 1^ 

where the subscript "la" has following convention: "1" stands for the first type subprocess, 

shown in Fig. Q "a" corresponds to diagram in Fig. ^a. In Eq. ([I9|), the factor Ci is an 

overall color factor for first type subprocess. The function Hia in Eq. (|I9[) is given by 

Hia = l^^ Tr [7 ■ (xV + krh " P7 ■ {x'p' + kT)Rtl -pLZ] {-9ap) {-g,.) , (20) 

where Ri^ and L"^ are the right and left blob, respectively, as shown in Fig. |^a. These blobs 
include all possible tree Feynman diagrams with the external partons shown in the figure. 
Substituting Eq. (|TUp into Eq. (|r^), we obtain 

H^^ = —Ci ^ I dxi dx2 dxe'^'P^'y e^^^''^'^ e^^^f^^^i"-^'^) 

An x's + u — J 

1 1 

X — 



^2 ^2 
^1 ~ + ^1 ~ ^2 - ^ - 

r f i - kt- 2kT ■ q\ - , , s 

x6\x + Xi-\ ■ ——] Hia{x,Xi,X2) . (21) 

y x's + u — Q^ J 

After performing dx by the 5-function, and dxi and dx2 by contour, we derive 

x's + u — 

x^(-l/2") Oiy" - y^) Hia{x, xi, X2) , (22) 

where the 6'-functions result from the contour integrations, and the momentum fractions for 
the function Hia are defined as 

i-2kT-q-kl k^ 

^ = ~, 7T2 T 23a 

x's + u — Q"^ x's 

xi = ^; (23b) 
x's 

X2 = . (23c) 
10 



Similarly, we derive contribution from the diagram shown in Fig. as 

^ ' x's + u-Q^ 

x^(2/^ - Vi) 0{y' - y^) Hib{x, xi, X2) , 

where x is defined as 



(24) 



X = —- 



(25) 



x's + u — Q'^ ^ 

and Xi and X2 are defined in Eqs. (|23b|) and ( p3c| ). Similar to Eq. (|20|) , the partonic part 
-ffib is given by 

1 



4 



7 ■ {x'p')R^^-f -pLZ i-9ap) {-gt^v) 



(26) 



The diagram shown in Fig. ^ has following contribution 

^ ' x's + u-Q^ 

xO{yi - y^) 0{-yi) Hibix, xi, X2) . 



(27) 



In deriving Eq. (P7|), we used the fact that the partonic part Hic = Hib when xi and X2 are 
evaluated at the same values as listed in Eq. (P5|). 

Combining i/ia, iJi;, and Hi, (given in Eqs. (p2l), (pif) , and (p7|), respectively) together, 
we obtain the total contribution to H, defined in Eq. i\L7\) , from the first type diagrams 
shown in Fig. ^, 



Hi = Hia + Hih + Hic 
1 

= (7ra J Ci 



X 



x's + u — 

^i-.p+y' H,,^^,X,,X2) - e'^'^'-y- Hiaix,Xi,X2) 



(28) 



All momentum fractions in Eqs. (^) are evaluated at the values defined in Eqs. (^) and 
(P5|). In deriving Eq. (pSj), we have dropped a term proportional to 
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This is because of the phase exp[ixp'^y~] which effectively restricts y~ ~ l/{xp'^) 0. 
Physically, it means that all y integrations in such term are localized, and therefore, will not 
give any large nuclear size enhancement. 

By substituting Eq. (|28| ) into Eq. (|T^), we can obtain the lowest order double scattering 
contribution from first type diagrams shown in Fig. ^ One important step in getting the 
final result is taking the derivative with respect to as defined in Eq. (0). Comparing 
Eq. ( p8|) with Eq. (^), and observing that 



^^,p+y- jj^^^-^ _ ^^.r,+y- ^^^^^^ 



kT=0 



, 



(29) 



we found that the derivatives on the exponential exp[i (/c|,/x's) — y^)] do not con- 

tribute, and that we can therefore set exp[i {k^ / x' s) {y{ — 2/2^)] = 1 in Eq. (p^). Substi- 
tuting Eq. (pHD into Eq. (0), we obtain 



da 



(SH) 



dQ'^dq^dy 



1 92 



X 



2 dk^dk^ 



2x's x's + u — Q"^ 



Tq{x, A) Hia{x, Xi, X2) - Tg{x, A) Hibix, Xi, X2) 



(30) 



with X, Xi, X2 and x defined in Eqs. ( P3| ) and (P5|). In Eq. (|30|) , Tq(x,y4) is the two-quark- 
two- gluon matrix element, and defined as 0: 



T,(x,A) 



27r 



271 



.e{y--y-)ei-y,) 



x-{pA\F^'-iy2)Moh^My~)F'-''iyi)\pA) 



(31) 



C. Final results for soft-hard processes 

The derivatives with respect to kx in Eq. (|30| ) are straightforward. It is most convenient 
to re-express derivatives with respect to kj- in terms of derivatives with respect to x or x . 
After working out the derivatives, we obtain the parton-nucleus cross section for first type 
soft-hard process 
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dcTi 



(SH) 



dQ^dq^dy ^« 

X 



1 



1 



2x's x's + u — Q'^ 



-Tg{x,A)Hi{x] 



+ 



dx"^ \x 
d_ n 

dx \x 



2g2 



Tg{x,A)H,{x) 



{x's + U- g2)2 



x's{x's + u — Q"^) J 

where x is given in Eq. (|25|) , and a^^ is the parton level Born Drell-Yan cross section 



(32) 



0", 



DY 



2 

em 



9Q2 



(33) 



In Eq. (|3^) , Hi{x) is the partonic hard part 



ifl(x) = CiXi^ia(x,Xi = 0,X2 = 0), (34) 

and ffia(x, Xi = 0, X2 = 0) is defined in Eq. ( pOD with fc^ = 0. In Eq. (^), we already 
include the factor (l/27r)(e^/3Q^) from the leptonic part of phase space (see Eq. ([TT|)). 

After convoluting Eqs. (|32D with the corresponding parton distribution from the beam, 
we obtain double scattering contribution to hA — > II from the first type soft-hard double 
scattering process 

dai^^^ sr^ f , , r / i\ d&i^^^ 



J2 J dx'fg/hix') 



(35) 



dQ^dqj^dy ^ J dQ'^dqj^dy ' 

where the parton level double scattering part d&i^^^ /dQ'^dq^dy is given in Eq. (p^). 

Following the same derivations, we obtain contributions from second type and third type 
processes shown in Fig. |^ and Fig. |^. For second type diagrams, the final result for second 
type soft-hard process contributing to hA II is 



dQHq^dy ^'^"^ ' dQHq^dy 



(36) 



where the parton level double scattering part da^^'' / dQ'^dq^dy is: 



da^^"^ 1 1 

dQHqldy " ^ 2x's x's + u-Q'' 



X ■ 



+ 



^^{^-Ux,A)H.ix) 
l(^T,ix.A)H.ix) 



2ql 



[x's + U- Q2)2 



u 



x's{x's + u ~ Q"^ 



(37) 
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In Eq. (|37|), the four-gluon matrix element Tg is defined as 



Tg{x,A) 



dy 

271 

1 



dyi dy2 

271 



.9{y~ - y-)9{-y,) 



{pA\F^+{y,)F^+{0)F^-{y~)F+{y^)\p^) 



(38) 



The partonic hard part H2 in Eq. (^) is defined as 



H2 = C2 H2a{x, Xi = 0, X2 = 0) 



(39) 



where C2 is the overall color factor for second type diagrams. H2a is given by the diagram 
shown in Fig. ^, and defined as 

1 



^^2a = -Tr 
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7 • {x'p + Xip + kT)R2a 7 ■ PZ LZ (-gap) (-^M^) 



(40) 



where = x'p' + {x + Xi)p — q is the momentum carried by the quark going to final state. 

Similarly, for the third type soft-hard process, as sketched in Fig. |^, the final result for 
third type soft-hard process contributing to hA — > II is 



da 



(SH) 



J2 J dx'fg/hix') 



da' 



(SH) 



dQ'^dqj'dy dQ'^dq^^dy 
where the parton level double scattering part da^^^^ / dQ'^dq^dy is: 



(41) 



2x's x's + u — Q"^ 



-Tg{x,A)Hs{x] 



+ 



dx"^ \x 
l{lT,ix,A)H,ix: 



2ql 



2(Q^-^) 
x's{x's + u — Q"^) 



(42) 



In Eq. (|4^), the partonic hard part i^s is defined as 



H2. = C3X H3a{x, Xi = 0, X2 = 0) , 



(43) 



where C3 is the overall color factor for the third type process, and H^a is given by the 
diagram shown in Fig. 0, and defined as 



(44) 
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where ps is the same as that defined following Eq. (^Ol). 

The partonic short-distance hard parts, defined in Eqs. (|3^), ( pOf ) and (|i3|), can be easily 
evaluated by calculating the corresponding Feynman diagrams. Our results are 

Hi{x) = I- (^^a, + -^a, +^^^^3) (45a) 

27 \ X X — Xa [X — Xa) J 



H2{x) = 4 (^3 + ^4 - ^^^^ 91 a,) (45b) 



-f^3(a;) = -7: f — — 0-5 + — ^5-6 , 0-4) , (45c) 



and 

= + (46b) 

-u 

= x's + u- ^^''^ 

a, = ^im^ (46d) 
x's 

^5 = (46e) 
cTg = — (46f ) 

In Eq. (|45|) , Xa is defined as 

In a special case, = 0, our result for soft-hard process corresponds to direct photon 
production [1^. In Eqs. (|3^), ([37| ) and (^^ , the only unknown quantities are the twist-4 
parton distributions defined in Eqs. (|31| ) and (|38D. They were first introduced by the authors 
of Ref. 0. These twist-4 parton distributions are universal functions, and as fundamental 
as the usual twist-2 parton distributions. They can be measured in one process and apphed 
to another process. More discussions about these twist-4 distributions will be given in Sec. 
VI. 
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IV. DOUBLE-HARD SCATTERING CONTRIBUTION 



In this section, we derive the contribution from double-hard processes. The calculations 
for double-hard processes are simpler than soft-hard processes. In soft-hard double scattering 
case, we use covariant gauge. In order to extract the contribution beyond leading-twist in 
the covariant gauge, we have to keep kx for the soft gluon in coUinear expansion and expand 
the hard-part at /ct = 0. For double-hard processes, since the leading behavior is twist-four 
contribution, we do not need to do such fcy expansion. In addition, only the symmetric 
diagrams contribute to double-hard scattering. At lowest order, there are three types of 
diagrams contributing to double-hard subprocesses, as shown in Fig. 0, Fig. § and Fig. |^. In 
the following subsections, we give the derivations for one subprocess and provide the results 
for other subprocesses. 



A. Separation of matrix elements and partonic hard part 

Consider subprocess shown in Fig. |^. Similar to soft-hard case, there are three indepen- 
dent four-momenta linking the partonic part and the corresponding two-quark-two-gluon 
matrix element. After collinear expansion, the parton- nucleus double scattering contribu- 
tion for subprocess shown in Fig. |^ can be written as: 



dQ'^dq^dy 2x's 



J dxdxidx2T^{x,xi,X2,p) H^^ {x'p',x,xi,X2,p,q) ■ (48) 



The superscript ^^DH" denotes double-hard scattering. The double-hard two-quark-two- 
gluon matrix element, T^, is defined as: 



t''{x,x„X2,p) = /^^^e--^^-e--^^re-«(---)^-^.- 
J 2tt 27r 271 
1 1 



y<--{PA\A"{y2MU^)i^Uy~)A'{y^)\PA) (-^i). (49) 

z p 

i~9ap) represents physically polarized gluons. In Eq. (|i8|), the partonic part H^^ is given 
by the diagram shown in Fig. |^ with gluon lines contracted with {—gP")/2, quark lines from 
the target traced with (7 ■ p)/2, and quark lines from the beam traced with (7 ■ {x'p'))/2. 
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In order to perform the integration of momentum fractions, we rewrite the double scat- 
tering contribution defined in Eq. as 



x' 
1 



dQdqj'dy 2x's J 2-n 27r 27r 



X - i,pAW{y-^. Ot) 7+ ^n{y-) A'ivi) \va) {-gtp)- (50) 

In Eq. ([50|), the modified partonic part H^^ is defined as 

xH'^" {x'p',x,xi,X2,p,q) , (51) 

All the integrals of the momentum fractions can now be done explicitly without knowing 
the details of the matrix elements. 

For diagram shown in Fig. 0, the final state photon-gluon two particle phase space can 
be written as 



r(2) = ^ 1 — s(x + x, + — 

IGvr^ x's + u — \ x'i 



t 



(52) 



'S + U- Q'^y 

Again, we neglect dQ^dq^dy in the F*^^^ because of the definition of the differential cross 
section. Using Eq. (|5^) , the contribution to H from the diagram shown in Fig. ^ can be 

expressed as 

-rjDH 1 ^ J i \ 



IGvr^ x's + M — \ x's + u — Q"^) 
X (^) 77^-, — H?''{x, xi, X,) . (53) 

\6/ {Q — U) X — Xa — %t X + X2 — Xa + 

where the subscript "1" stands for first type double hard scattering subprocess shown in 
Fig. 1^. Xa is defined in Eq. (^). The function H^^ in Eq. (|53D can be represented by the 
diagram shown in Fig. ^a. It is given by 

^ = ^ Tr [7 ■ xy/^f 7 ■ {q - xip) L7] {-g^p) {-g^,) , (54) 

where R^" and L^^ are the right and left blob, respectively, as shown in Fig. |^ These 
blobs include all possible tree Feynman diagrams with the external partons shown in the 
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figure. The factor 1/3 in Eq. (|53|) comes from averaging over the color of the quark which 
annihilates to produce 7*. We write it out explicitly here so that the remaining color factor 
will be just the color factor for H^'^ , which is represented by the diagram in Fig. llOa. 



Substituting Eq. (|5^ ) into Eq. (0), we perform dxi integration by using the 5-function, 
and dx, dx2 by contour integration. We then have 

fjDH _ ]_ ^ ixap+y~ Jxip+{y~-y^) 

^ 12 x's + u- g2 (g2 _ ^) 

xei-y^)9iy--y^)H^^{x,x,,X2) , (55) 

where the 6'-functions result from the contour integrations, and the momentum fractions for 
the function are 

x = Xa = -p^ — (56a) 

Xi=Xb= — — — - Xa (56b) 

x's + u — Q'^ 

X2 = (56c) 
Substituting Eq. (pS]) into Eq. (0), we obtain 



d ^fA% ^ ( 87r2a|a^\ J_ r dy_ {p+fdy^ dy^ ixaP+y-ix,p+(y--y-) 



dQ^dq^dy \ 3 J 2x's J 27i 2n 

X I j^^, {pA I (y^- , Ot) (0) 7+ V^, (y- ) F^'^ (yr ) \pa) 

xei-y^) e{y- - yf) , J — ^i^^^(a;, Xi, x^) . (57) 
x's + u — Q"^ [Q"^ — u) 

In Eq. (|57|) , we already convert the field A'^ into field strength F^'^ by partial integration. 
Define the double hard matrix element [01 



TS,{xa.x,,A) = ^/^/^/ P^dy:^e{y- -y-,)e{-y,)e^^^^'-y- e^^^^^y^-y-^^ 



J_ f dyi f dy- 

XbP^ 

^{PAK^{y2)Uoh^My')F^''iyi)\pA) ■ (58) 



Substituting Eq. (|58D into Eq. (^Tj), we obtain the parton level cross section for first type 
double hard process 
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1 _ 2 DY / A 2\ 

dQ^dq^dy ~ ^ """^''x's + u-Q^ {Q^- u) 

x^^.?(^-^^'^)4''''' (59) 

with cr^^ the parton level Drell-Yan Born cross section defined in Eq. (p3|), and partonic 
hard part H^^ given in Eq. (|5^). In deriving Eq. (|59D, we already included the factor 
(l/27r)(e^/3(5^) from the leptonic part of phase space. 

Comparing Eq. ( |59D with Eq. (|32D, we see that the contribution from double-hard process 
does not have the derivative terms, while the contribution from soft-hard process depends on 
the derivatives of the distributions and hard part. The derivatives of x in Eq. ( P^ ) actually 
result from the derivatives with respect to /ct in Eq. (|30D . As we pointed out in Sec. Ill- A, for 
soft-hard process, the leading term when kx = is the leading twist eikonal contribution and 
does not correspond to physical double scattering. The lowest order high-twist contribution 
comes from the second order derivative term in expansion. For double-hard process, 
since momentum fraction of both partons are finite, the leading contribution is already a 
twist-four contribution. 



B. Final results for double-hard processes 



Convoluting Eq. ( [59| ) with the corresponding parton distribution from the beam hadron, 
we obtain contribution to hA — * // from first type double-hard process shown in Fig.|^: 

dQVrdy ^ ? i '^'"'^'"'^'''^Wd^ ' 

with dai^^^ /dQ'^dq^dy given in Eq. (|^). 

Following similar derivations, we obtain contributions from second and third type double- 
hard processes shown in Fig. |^ and Fig. |^. For contribution from second type double-hard 
process shown in Fig. ^, we have 

dQUq\dy '^'"'^'"'^'''^dQk^d^ ■ ^^^^ 
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The parton level double scattering part da2^^^ / dQ'^dq^dy is given by 



da. 



(DH) 



-.2 ^DY ( 2\ 

do (47raJ — 



dQ'^dqj'dy 



1 



1 



X'S + U- Q2 (Q2 _ ^) 



n I T iXa,Xb,A)H2 

2x'XbS 



(62) 



cr(f^ is defined by Eq. (^). In Eq. (|62D, the partonic hard part H2^^^ is given by 



H. 



(DH) 



4^ 



7 ■ psR^"" l-{q- XaP) L^^] {-9al3) {-9ap) 



The diagram representing H2^^^ is shown in Fig. 10b. 

For the third type double-hard process, shown in Fig. we have 



da. 



(DH) 



dQ'^dqj.dy 



J dx'fg/h{x'] 



da. 



(DH) 



dQ'^dqj^dy 



And da^^^ /dQ'^dq^dy is given by 



dQHq^dy 9 ° ^ x's + u-Q^ {Q^ -u) 



ZX'XbS 

In Eq. (|65|) , the partonic hard part H^^^^ is given by 



(DH) 



(63) 



(64) 



(65) 



H. 



(DH) 



-Tr 



(66) 



The diagram representing H^^^^ is shown in Fig. [I0| c. In Eq. (|65|) , the four-quark matrix 
element is defined as 

T.'iix^^Xb^A) =J^J^J P^dy^e{y--y^)e{-y,)e^^^^'-y-e^^^^'-'~y^'y^^ 

x{pA\My2h^Myi)M^h^My')\PA) ■ (67) 

The partonic short-distance hard parts h[^^\ h'^^^ and h'^^\ defined in Eqs. (0), 
and (§^), can be evaluated by calculating corresponding Feynman diagrams in Fig. [1^. 
The results are: 
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friDH) _ 4 f x's Q^-U \ {X'sf + {Q^ - uf 

"el^ Q^-u ^ x's + u-Q^J 8 (a:'s)2 ' 

Eq. ( |68|) already includes the color factor for each individual Feynman diagram. 

In above derivations, the only unknown quantities are double-hard matrix elements T^^ 
and T^, defined in Eqs. (|58D and (67). These double-hard matrix elements are universal 



functions. They cannot be calculated perturbatively, but can be measured in one process 
and tested in another process. Theoretical predictions for double scattering can also be 
tested once these matrix elements are measured. 

V. RELATION BETWEEN SOFT-HARD AND DOUBLE HARD SCATTERING 

As discussed in Sec. II, the contributions from soft-hard processes and double-hard pro- 
cesses correspond to the contributions from different residues when we integrate over the 
parton momenta using contour integration. Let us now look at the relationship between 
the two kinds of double scattering process. As an example, we exam the process shown in 
Fig. ^. The full amplitude for the double scattering has the general form: 

r I 1 

M ^ dxi F(4,a;i,x), (69) 

J Xi — Xia + ie Xi — Xib + ie 

where Xi is the parton momentum fraction which needs to be integrated, and x is the other 
parton momentum fraction. The sum of these two parton momentum fractions is fixed by 
the kinematics, such as x'p', and q. In Eq. (|69D, the pole xi —xia corresponds to the first, 
and Xi — Xib corresponds to the second propagator, as shown in the figure. The poles are 
given by 

xia = — , (70a) 

+ 2p3 ■ kr + 2p3 ■ x'p' 

xib = ; ^ , 70b 

x's — 2p3 ■ p 
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and they are both in the same half of the complex plane in the region of interests {x's — 2p3 ■ 
p > 0). In Eq. (|69D, the function F{k^,Xi,x) is a non-vanishing and smooth function when 
Xi = Xia and xi = xu- It is proportional to the parton fields of momentum at xip and xp. 

When Qt ^ and — > 0, xia —>■ while xib is finite. Using contour integration to carry 
out dxi-integral in Eq. (p5|), we have 

F{k^, Xia, Xtot — Xia) _ -^(^T; Xib, Xtot ~ Xib) 
X\a X\}y X\a X\}) 

—hard —hard j 

(71) 



M 



where Xtot is the sum of the total mometum fraction from the target, and is a function of x', 
Ps and q. Note that in Eq. (|7lD, the amplitude of these two terms have the opposite sign. 
The square of these two terms gives the leading high-twist contribution of double scattering. 
Square of the first term corresponds to the soft-hard double scattering case, and square of 
the second term corresponds to the double-hard double scattering case as we derived in the 
last two sections. As we will discuss later, the interference term is small when qj- is large. 
When qx = 0, p^^ = kx, and the two poles given in Eq. ([70D become 



1^2 

= ^ , (72a) 

xi, = x-v/i/(2p3.) - 1 ^2 ^ _ (72b) 

x's - 2A/ip32 - ^/s k^/ (2p3^) ^ 2^p3^ 

As /ct — ^ 0, both xia and xu vanish. To understand the region where ~ and fcy ~ 0, all 
terms proportional to Xi (or xia, or xif,) should not be included in F{kx,xi,x = Xtot — xi) 
for the leading pole contribution. Therefore, for the leading pole contribution in the region 
of ~ and /c^ ~ 0, F ~ F{xtot) and 

—hard —hard ■ 

(73) 

X\a XH) 

Consequently, M ~ Mgoft-hard — M^oubie-hard gives no leading contribution to double scat- 



tering when qT = [11 



The double-hard process corresponds to the classical double scattering picture. The 
contributions of double-hard processes given in Eqs. (^Of), (|62D and (|65| ) have the form of 

22 



c7 ~ (7i • (^2 ■ , (74) 

where ai and (J2 are the partonic cross section for first and second scatterings respectively. 
The soft-hard process does not have the corresponding classical double scattering picture, 
and hence the cross section does not have the form of Eq. (|7^) . The fact that the scattering 
amplitude for double-hard process and soft-hard process cancel each other when = is a 
pure quantum effect. 

From Eq. (^), we can see that the interference term of soft-hard double scattering 
and the double-hard double scattering will have a negative sign. That will give "suppres- 
sion" effect. The interference is proportional to F*{k^,xia,Xtot — Xia)F{k^, xib, Xtot — xu) 
or its complex conjugate. As kr 0, F*{k^,xia,Xtot — xia)F{k^,xib,Xtot — xu) —>■ 
F*{0,0,xtot)F{0,xib,Xtot — xib). If xib 7^ 0, overlap in phase space for F* and F is clearly 
small because of the difference in parton momenta, while full overlap takes place when 
Xib 0. Therefore, in large Qt region, we only consider the contribution from soft-hard 
and double-hard process respectively, and neglect the interference between them. However, 
in smaller qt region, one of the two hard partons in the double-hard amplitude will become 
softer, and may have much bigger phase space to overlap with the soft-hard amplitude. 
Consequently, the interference term may become important and give large suppression ef- 
fect. In order to obtain correct qt dependence for smaller qt region, we need to consider the 
interference between soft-hard and double-hard scattering as well as doing resummation of 
powers oi asln'^iQ'^/qT) 0- 

VI. NUMERICAL RESULTS AND DISCUSSIONS 

In this section, we present our numerical results for the nuclear enhancement in Drell- 
Yan transverse momentum distribution. We numerically evaluate the ratio of differential 
cross section to single scattering contribution R defined in Eq. by using our analytical 
results presented in Sec. Ill and Sec. IV. 
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The ratio R defined in Eq. (|]) depends on contributions from botli single scattering and 
double scattering. All these contributions depend on the non-perturbative parton distribu- 
tions or multi-parton correlation functions. In deriving following numerical results, the Set-1 



pion distributions of Ref. are used for pion beams; and the CTEQ3L parton distribu- 
tions of Ref. |TB| are used for free nucleons. The twist-4 multi-parton correlation functions 
defined in Eqs. (0), (|38D , (|58|) and (|^) have not been well-measured yet. By comparing the 



definitions of the soft-hard correlation functions Tq and Tg with the normal twist-2 parton 
distributions Jill, authors of Ref. [H proposed following approximate expressions for the 



twist-4 soft-hard correlation functions Tq and Tg-. 

Ux,A) = X^A^I^h,A{x,A) (75) 

where i = q,q, and g. The fi/A are the effective twist-2 parton distributions in nuclei, and 
the factor A^^^ is proportional to the size of nucleus. The constant has dimensions of 
[energy]^ due to the dimension difference between twist-4 and twist-2 matrix elements. The 
value of should be determined by experimental measurement. It was estimated in Ref. 0] 
by using the measured nuclear enhancement of the momentum imbalance of two jets in 
photon-nucleus collisions |]IR|TB|, and was found to be order of 



A^ ~ 0.05-0. IGeV^ . (76) 

This value is not too far away from the naive expectation from the dimensional analysis, A^ ~ 
^QCD- The theoretical predictions for the nuclear dependence in direct photon production 
base on this value was consistent with the experiment |]10|- In our calculation below, we use 
A2 = 0.05 GeV^. 

From the definition of the correlation functions in Eqs. ( pT)) and (|38|), the lack of os- 
cillation factors for both yjf and y2 integrals can in principle give nuclear enhancement 
proportional to A^^^. The A^^^ dependence is a result of the assumption that the positions 
of two field strengths (at and y^, respectively) are confined within one nucleon. 

For double-hard matrix elements Tgg and Tgg defined in Eqs. ( pBD and (0)) adopt the 



model of Ref. and take 
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T,g{xa,xi„A) = ^A^r^f,{xa)Ux,) ' (77a) 
T,,(x,,a:,, A) = ^A'l^f,{xa)h{x,) , (77b) 

where fg, fq are normal twist-2 quark and anti-quark distribution, and fg is a corresponding 
gluon distribution. C = (0.35/87rrQ) GeV^. vq is the value of nucleon radius, and tq ^ 
1.1 — 1.25. Here we already convert the unit for C to GeV^. The factor (l/27r) comes from 
the difference in the overall factor of the definition for matrix elements in this paper and in 
Ref. p. 



In Eq. (|75D, the effective nuclear parton distributions fi/A should have the same operator 
definitions of the normal parton distributions with free nucleon states replaced by the nuclear 
states. For a nucleus with Z protons and atomic number A, we define 

f,/Aix, A) = A (^/./^(x) + |/./p(x)) i?™^(x. A) , (78) 

where fi/Nix) and fi/p{x) with i = q,q,g are normal parton distributions in a free neutron 
and proton, respectively; and N = A — Z . The factor i?™^ takes care of the EMC effect 
in these effective nuclear parton distributions. We adopted the R^^'^ from Ref. [T^, which 
fits the data well. 



In Fig. |lT], we show the ratio of double scattering and single scattering versus qt- The 
ratio is normalized by A^^^. The overall feature is that the ratio becomes smaller for larger 
Q and qt- This is the typical behavior for high twist contribution. It is suppressed by the 
factor of l/Q"^ or Q and qt are both large scales in the processes. On the other 

hand, we see that when qjn becomes very large, the ratio starts to rise again. This is the 
effect from the edge of phase space. The soft-hard double scattering contribution depends 
on the derivatives of the matrix elements. At edge of the phase space, the derivative term 
become more important. When x ^ 1, the parton distribution f{x) has (1 — x)° behavior, 
with a > 1. As X — 1, df{x)/dx » f{x). Therefore, the double scattering contribution 
becomes more important at the edge of the phase space. The curve for Q = 11 GeV rise 
more steeply than the curve for Q = 5 GeV, because at high Q, we approaches edge of phase 

25 



space faster as we increase qr- 



Fig. |T2| shows the ratios of the soft-hard and the double-hard contributions to single 
scattering contribution. The ratio is normalized by A^^^. We see that the size of soft-hard 
and double-hard processes are comparable. At higher qx, the ratio for soft-hard process 
increase, as explained above. The contribution from double-hard processes increase faster 
when qT decreases. The contribution from double-hard processes depends on the product 



of two parton distributions as given in Eq. (|77b|) . The momentum fractions Xa and Xb are 



smaller than the momentum fraction for soft-hard processes. Because of the steep increase 
in parton distribution for smaller momentum fraction, the contribution from double-hard 
processes increase faster as qx becomes smaller. 



Fig. shows the comparison of double scattering contribution and single scattering con- 
tributions. The single scattering contribution is normalized by A and the double scattering 
contribution is normalized by A^^^. Without extra A^^^ factor, the double scattering contri- 
bution are small compare to the single scattering contribution. Perturbative calculation for 
double scattering should be reliable. Because of A^^^ enhancement, the double scattering 
becomes important. As xp increases, both double scattering contribution and single scatter- 
ing contribution decreases, but double scattering contribution decreases slower than single 
scattering. We expect larger enhancement in larger xp region. 



In Fig. |14|, we show the nuclear dependence of the qx spectrum in large qx region for 
different nuclear target: C, Ca, Fe, and W, when 800 GeV proton beam is used. The 
figure is plotted for R versus qx- The quantity R is defined in Eq. (|^). It is the ratio of 
total differential cross section and the single scattering cross section. Because of the A^^^ 
enhancement for double scattering contribution, the ratio is larger for larger A. 

As we discussed in the last section, when qx is smaller, the interference between the 
double hard scattering and the soft-hard scattering becomes important. The interference 
will give suppression. We expect that the double scattering contribution will decrease as 
qx decreases in smaller qx region. Fig. K^ shows our expected ratio of double scattering 



contribution to single scattering contribution. The dotted curve and the curve on the right 
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side of the dash hne is obtained by using our results presented in this paper. The curve 
on the left side of the dash line is our expected ratio when becomes smaller. Since 
the integrated cross section da/dQ"^ has little A-dependence, and we have shown nuclear 
enhancement in large qt region in this paper, we expect that there will be suppression in 
small qt- In small g-r region that the interference between the soft-hard and double-hard 
processes becomes dominant contribution. From the figure, we see that the calculated ratio 
starts to rise steeply when qt approaches 3 GeV and below. This is the sign where the large 
logarithms of agln'^iQ'^ /qt) and the interference become important. 

In summary, both soft-hard and double-hard double scattering processes contribute to 
the nuclear dependence in Drell-Yan pair production. The double-hard subprocess resembles 
the classical double scattering picture, while the soft-hard subprocess and the cancellation 
of the leading double-hard and soft-hard contribution are pure quantum effect. Both types 
of double scattering processes give A^^^ nuclear enhancement. However, the interference 
between these two types of process have the suppression effect. In large qx region, the inter- 
ference can be neglected and we have A^^^ nuclear enhancement for Drell-Yan qx spectrum. 
In smaller region, the interference term becomes more important and may give nuclear 
suppression. Calculation of nuclear dependence of Drell-Yan process in small qx region will 
be addressed in a future publication |p. 
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FIGURES 

FIG. 1. Diagram for single scattering: a) Annihilation diagram; b) "Compton" diagram. 

FIG. 2. A Graphical representation of double scattering contributions in proton- nucleus colli- 
sions. 

FIG. 3. Sample Feynman diagrams contributing to next-to-leading order double scattering: (a) 
soft-hard process; (b) double-hard process. The circles represent the poles when we integrate over 
the parton momentum. 

FIG. 4. First type soft-hard double scattering subprocesses: "Annihilation" diagrams corre- 
sponding to the two-quark-two-gluon matrix element; a) real diagram, b) and c) interference dia- 
gram. 

FIG. 5. Second type soft-hard double scattering subprocesses: the real "Compton" diagrams 
corresponding to the four-gluon matrix element. 

FIG. 6. Third type soft-hard double scattering subprocesses: the real "Compton" diagrams 
corresponding to the two-quark-two-gluon matrix element. 

FIG. 7. Diagram corresponding to first type double-hard subprocess. 

FIG. 8. Diagram corresponding to second type double-hard subprocess 

FIG. 9. Diagram corresponding to third type double-hard subprocess 

FIG. 10. Diagrams representing the partonic hard parts for double-hard processes 
for double-hard process shown in Fig. |^; b) 62'^ for double-hard process shown in Fig. [ 
for double-hard process shown in Fig. |^. 
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FIG. 11. The ratio of double scattering contribution and single scattering contribution normal- 
ized by A^^^ for a SOOGe^ proton beam when Q = 5GeV and llGeV. 

FIG. 12. The ratio of soft-hard and double-hard contributions to single scattering contribution. 
The ratio is normalized by A^^^. 

FIG. 13. Comparison of double scattering and single scattering contributions. 

FIG. 14. A-dependence of qt distribution at large qt region for nuclear targets C, Ca, Fe 
and W. The quantity R is the ratio of total differential cross section and the single scattering 
contribution, as defined in Eq. @. 

FIG. 15. Expected ratio of double scattering to single scattering when qx is small. 
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